MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
ARTIFICIAL INTELLIGENCE LABORATORY 


A.I. Memo No. 1332 November, 1991 

The Alignment of Objects with Smooth Surfaces: 

Error Analysis of the Curvature Method 

Ronen Basri 


Abstract 

The recognition of objects with smooth bounding surfaces from their contour images is 
considerably more complicated than that of objects with sharp edges, since in the former 
case the set of object points that generates the silhouette contours changes from one view 
to another. The “curvature method”, developed by Basri & Ullman [1988], provides a 
method to approximate the appearance of such objects from different viewpoints. In this 
paper we analyze the curvature method. We apply the method to ellipsoidal objects and 
compute analytically the error obtained for different rotations of the objects. The error 
depends on the exact shape of the ellipsoid (namely, the relative lengths of its axes), 
and it increases as the ellipsoid becomes “deep” (elongated in the Z-direction). We show 
that the errors are usually small, and that, in general, a small number of models is 
required to predict the appearance of an ellipsoid from all possible views. Finally, we 
show experimentally that the curvature method applies as well to objects with hyperbolic 
surface patches. 
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1 Introduction 


Visual object recognition requires the identification of objects observed from different 
viewpoints. In the alignment approach [Ullman 1989] objects are recognized in a two 
stage process. First, the pose of the object is determined, and then the appearance of 
the object from this pose is predicted and compared to the image. Such a method is 
used in [Fischler Sz Bolles 1981, Lowe 1985, Faugeras Sz Hebert 1986, Chien Sz Aggarwal 
1987, Huttenlocher &: Ullman 1987, Lamdan et al 1987, Thompson & Mundy 1987]. In 
these works, the method is applied to either planar or polyhedral objects. 

Basri & Ullman [1988] extended the alignment approach to handle rigid objects 
bounded by smooth surfaces. The method is called “the curvature method”, and it 
approximates the appearance of such objects from different viewpoints using the 3-D 
curvature of points along the contours. The authors showed experimentally that in gen¬ 
eral a few models are sufficient to predict the appearance of objects from all possible views 
with high accuracy. This approximation is also the key for showing that the appearance 
of objects with smooth bounding surfaces can be predicted by linearly combining a small 
number of views [Ullman Sz Basri 1991]. 

In this paper we analyze the curvature method. We apply the method to ellipsoidal 
objects and derive an expression that describes the errors obtained. We analyze this 
expression and show that the error depends on the exact shape of the ellipsoid (namely, 
the relative length of its axes). We show that the errors are usually small, and that, in 
general, a small number of models is required to predict the appearance of an ellipsoid 
from all possible views. Finally, we show experimentally that the curvature method is 
not restricted to contours generated by parabolic surface patches, but it can also handle 
contours generated by hyperbolic patches. 


2 The Curvature Method 

In this section we briefly review the curvature method. The difficulty in predicting the 
appearance of objects with smooth bounding surfaces is described below. The silhouette 
of an object is the set of its bounding contours observed in the image. These points 
are generated by the rim. of the object [Koenderink & Van Doom 1981] (also called the 
contour generator [Marr 1977]), which is the set of all the points on the object’s surface 
whose normal is perpendicular to the visual axis. When the object contains only sharp 
edges, these edges compose the rim in all views of the object. The silhouette then includes 
the projection of those edges that are visible. When, however, the object is bounded by 
smooth surfaces, the rim is not fixed. Instead, it smoothly changes its position on the 
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object with viewpoint. To illustrate this phenomenon imagine for instance the shape of 
an egg. When the observer moves slightly, the observed silhouette of the egg is generated 
by a new set of points along the egg shell, and thus the silhouette appears to change its 
shape. 

To accurately align objects with smooth surfaces one must account for this special 
property of the objects. Previous approaches to this problem included the construction 
of detailed 3-D descriptions for the objects using either simple volumetric or surface 
primitives [Marr Sz Nishihara. 1978, Brown 1981, Dane Sz Bajcsy 1982, Potmesil 1983, 
Brady et al 1985, Faugeras Sz Hebert 1986], or using 3-D wires [Baker 1977]. These 
approaches, when handling objects with relatively complex structures, often require large 
memory to represent the objects and perform extensive computations to predict their 
appearances. 

Basri Sz Ullman [1988] suggested a novel approach to predict the appearance of objects 
with smooth surfaces. In their scheme an object is represented by its silhouette as seen 
from a particular viewpoint. Every point along the silhouette is associated with its spatial 
coordinates together with its 3-D curvature. The appearance of the object from other 
viewpoints is then predicted simply by applying linear transformations to the object. 
The authors showed experimentally that this method accurately predicts the appearance 
of objects for relatively large transformations. 

The curvature method is based on the following observation. The position change 
of the rim on the object depends mainly on the 3-D curvature of the object at the rim 
points. When the curvature is high the rim changes only slightly if at all. (In case the 
rim lies on a sharp edge, the curvature is infinite and the rim is fixed.) The lower the 
curvature is, the more significant is the position change of the rim on the object. 

Following this observation the curvature can be used to approximate the position 
change of the rim. The basic idea is illustrated in Figure 1. Consider an object 0 rotating 
by a rotation R around the vertical axis, Y. Let p be a point on the object’s rim. The 
figure shows a horizontal section of the object through p. Let r x be the curvature radius 
of p in this section, and let r x be a vector of length r x parallel to the X-axis. When the 
object rotates by R, point p ceases to be a rim point, and it is replaced by a new point, 
p ', approximated by 

p « R(p — r x ) + r x (1) 

This approximation holds as long as the circle of curvature provides a good approximation 
to the section at p. The method handles general 3-D rotation simply by substituting the 
vector r x in Eq. (1) with the radius vector that corresponds to the rotation applied. The 
radius vectors of a point with respect to all 3-D rotations can be recovered from a single 
number, the magnitude || (r x , r„) ||, where r x and r y are the curvature radii with respect 
to the vertical and the horizontal axes respectively. This number, together with the 
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Figure 1: The curvature method, (a) A horizontal section of an ellipsoid, p is a point on the rim, r is 
the radius of curvature at p, o is the center of the curvature circle, and a is the intersection of the Y-axis 
with this section, (b) The section rotated, p' is the new rim point, and it is approximated by eq. (1). 
(Borrowed from [Basri h Ullman 1988]). 


spatial coordinates of the point, is the only information required to apply the curvature 
method. 

Consequently, the authors suggested to represent an object by its contour image, 
as observed from a particular viewing direction. Each point along the silhouette has 
associated with it, along with its spatial coordinates, the magnitude of its curvature 
vector. To apply a rotation to the object, the appropriate curvature radius should be 
computed, and the new position of the point in the image is determined using Eq. (1). 
Translation and scaling are applied to the object in a straightforward manner. The 
method was implemented and applied to a number of objects, and it was shown that a 
few models can predict the appearance of these objects from all possible viewpoints with 
high accuracy. (See an example in Figure 2.) In addition, Ullman h Basri [1991] showed 
that the curvature information is implicit in a small number of views of the object. They 
used this observation to predict the appearance of objects with smooth bounding surfaces 
by linearly combining a small number of images. 


3 Properties of the Curvature Method 

The appearance of objects with sharp boundaries (for which the radius of curvature 
is zero) and of spherical and cylinderal objects is predicted exactly by the curvature 
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(a) (b) (c) 

Figure 2: Application of the curvature method to VW cars, (a) A contour image of a VW car. (b) A 
prediction of the appearance of the car. (c) Matching the prediction to the actual image (an overlay of 
(a) and (b)). (Borrowed from [Basri & Ullman 1988]). 

method. The appearance of smooth objects with arbitrary structures is, however, only 
approximated by this method. In order to demonstrate the properties of the curvature 
method, we applied this method to ellipsoids and analyzed the errors obtained. The 
analysis is given in this section. We first compute the errors obtained when a canonical 
ellipsoid rotates around the vertical (F) axis. We then compute the errors obtained when 
the same ellipsoid rotates arbitrarily in 3-D space, and show that the errors obtained in 
the two cases are similar. The error depends on the shape of the ellipsoid, in other 
words, on the relative length of its axes, and it increases as the ellipsoid becomes “deep” 
(elongated in the Z-direction). We show that the errors are usually small, and that, in 
general, a small number of models is required to predict the appearance of an ellipsoid 
from all possible views. 

We start with a brief explanation of the error function used. Consider an ellipsoid 
rotating about some axis V in the image plane. Let pi = (xi,yi) be the projected location 
of some rim point. Following rotation, the rim changes, and the point p x is replaced by 
a new point, p 2 = (a: 2 , 2 / 2)5 such that the vector p 2 — P\ is perpendicular to V . Denote 
the approximated location of p 2 according to the curvature method by p 2 = (x 2 , y 2 ). The 
observed error is measured by || p 2 — p 2 ||. Clearly, if we scale the ellipsoid the observed 
error would scale as well. We therefore need instead to consider a relative value for the 
error that is independent of scale. 

We define the error as follows. Consider the planar section through p 1 that is per¬ 
pendicular to the rotation axis V. This section forms an ellipse (or a single point in case 
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of a tangential plane). Let p 0 = (x 0 ,y 0 ) be the center of this ellipse. The relative error 
is defined by 

F II h - P2II 

II Pi ~ PO II 

E reflects the observed error relative to the projected size of the ellipsoid. Notice that E 
is independent of translation and scale of the ellipsoid. 


3.1 Rotation Around the Vertical Axis 


Let 


x 

— T 

a 1 


y_ 

b 2 


+ = 1 


be the surface of a canonical ellipsoid. Let p\ = ( X\,y \) be a point on its silhouette. 
When the ellipsoid rotates about the vertical (F) axis by an angle 0, pi disappears and 
is replaced by a new contour point p 2 = (x 2 ,y 2 ) with an identical y- value, y 2 = y\. Let 
p 2 = (x 2 ,y 2 ) be the approximated position of p 2 according to the curvature method. The 
horizontal section of the ellipsoid through p\ is an ellipse centered around p 0 = (0 ,y 0 ). 
Notice that the points p 1? p 2 , P 2 , and p 0 all lie on the same horizontal section implying 
that yi z= y 2 = y 2 = y 0 . The relative error is therefore reduced to 


E 


x 2 - x 2 


X\ 


(For reasons of convenience we ignore the absolute value operation in the discussion 
below.) 

Proposition 1: The error is given by 

(? c? I ? 

E( —, 6) = cos 0 H —-(1 — cos 0) — \ cos 2 0 -j—- sin 2 9 
a z a z V a 2 

(A proof is given in Appendix A.) 

The expression obtained for the error depends on two parameters the aspect ratio of 
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the ellipsoid, ^ 7 , and the angle of rotation, 9. Consequently, the error is invariant under 
a uniform scaling of the ellipsoid. 



3.2 Properties of the Error 

The prediction error obtained by the curvature method for a canonical ellipsoid rotating 
around the F-axis vanishes in the following three cases 
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Error Error 



Figure 3: The errors of the curvature method as a function of 9, the angle of rotation, (a) = 91 

\ and jp (b) £ = 2, 4, 9 and 16. (The parameters correspond the curves at increasing heights.) 


• 9 = 0 (that is, no rotation). 

2 

• ^ 2=1 (that is, c = a, the cross section is a circle). 

2 

• J 2 = 0 (that is, c = 0, a planar ellipsoid). 


As a function of 0, the angle of rotation, the error function is symmetric, that is, 
similar errors are obtained both for positive and negative angles. The absolute value of 
the error increases monotonically with the absolute value of 9. The partial derivative Eg 
also changes monotically with 9 , so the error increases faster for larger values of 9. The 
derivative is given by 

^ / c 2 . . „. cos 9 

Eq = (1--) sin 9(—j= ..■ = - 1) 

a \J cos 2 9 -f ^7 sin 2 9 



and assumes the following values 


• E e ( 0°) = 0. 

• E e ( 90°) = ^-1. 

Figure 3 shows the error as a function of 9 for several ellipses. 
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Figure 4: The maximal value of the error for canonical ellipsoids with c < a as a function of 9, the 
angle of rotation. 

As a function of ^ 7 , the relative size of the axes of the ellipsoid, the error behaves 
differently in each of the two ranges: (1) when c < a, and (2) when c > a. In the first 
case the ellipsoid’s width is larger than its depth. The error assumes small values even 
for fairly large values of 9. The maximal error is obtained when 

_ 3 1 

a 2 4 2(1 -f cos 9) 

and it assumes the following values 

• 0.24% at 30° (£ = 0.482). 

• 1.26% at 45° (£ = 0.457). 

• 4.14% at 60° 4 = 0.417). 

Figure 4 shows the maximal error as a function of 9. 

When c > a, the ellipsoid is deeper than it is wide, the error assumes larger values and 
is unbounded when 9 increases to 90°. The partial derivative E c 2 increases monotonically 

2 . 2 

with ^7 and reaches its maximum when ^7 —» 00 , where the error increases linearly in 
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Figure 5: The error as a function of the aspect ratio, of the ellipsoid for 6 = 20°, 30°, 45°, and 
60°. (The parameters correspond the curves at increasing magnitude.) 


%. The derivative is given by 


E c 2 = (1 — cos 0) 


7? 


and assumes the following values 


sin 


2 0 


2^/cos 2 9 + sin 2 0 


• £ e *(0) = - 




2 cos 0 


£ c 2 ( 1 ) = 




• lim c 2 E c 2 = 1 — cos0. 

—n* —► OO —TT 

rt ^ n* 



A model for such an ellipsoid would therefore cover only a restricted range of rotations. 

Larger rotations should be treated by additional models. Figure 5 shows the error as a 

2 

function of % for several values of 0. 

When a complete set of models is prepared for the appearance of an ellipsoid to be 

predictable from all possible views, it should be considered that following a rotation of 

90° about the F-axis, a and c, the axes lengths of the ellipsoid, interchange their roles. 

Therefore, an ellipsoid with c < a changes after a rotation of 90° to an ellipsoid with 

2 

c> a. An ellipsoid with a high aspect ratio, ^, changes to an ellipsoid with a low aspect 
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£i 

a 2 

1% 

2% 

3% 

4% 

5% 

6% 

2 

3 

3 

2 

2 

2 

2 

4 

3 

3 

3 

3 

2 

2 

9 

4 

3 

3 

3 

3 

2 

16 

4 

3 

3 

3 

3 

2 

49 

4 

3 

3 

3 

2 

2 

100 

4 

3 

3 

2 

2 

2 


Table 1: Number of models as a function of allowed error. 

ratio. Consequently, the small range of rotations covered by a model for an ellipsoid with 
a high aspect ratio is compensated by the large range of rotations covered by a model for 
the same ellipsoid after a 90° rotation. A small number of models is therefore required 
to represent the ellipsoid from all possible views. 

Table 1 shows the number of models required to cover the entire range of rotations 
about the T-axis for several ellipsoids. Because of symmetry considerations only rotations 
up to 90° should be considered. We see from the table that this number is small and 
does not exceed four even for extreme aspect ratios and allowed error of 1%. 

In preparing this table each ellipsoid was initially represented by two models, one 
taken at its canonical position, the other following a 90° rotation. If the two models did 
not cover the entire range of rotations, additional models were added at intermediate 
positions. In this case the value of the error is somewhat different then the canonical 
case. An expression describing this value is given in Appendix A. 


3.3 Rotation in 3-D Space 

We now consider the case of a canonical ellipsoid rotating arbitrarily in 3-D space. A 
rotation in 3-D space can be decomposed into three successive rotations, about the Z-, 
Y-, and Z-axes. The last rotation can be ignored since it does not deform the image and 
therefore does not change the errors. Let 

(—a: sin a -f y cos a 

b 2 

be the surface of a canonical ellipsoid rotated about the Z-axis by an angle a. We now 
examine this ellipsoid as it rotates about the T-axis by an angle 0. 



(x cos a + y sin a)' 


+ 


a* 
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Proposition 2: 


The error is given by 


C 2 

E <*••> 


in Eq. (2), where 


C 2 c 2 2 c 2 . 2 

= — cos a -f — sm a 


A 2 


a 1 


b 2 



(A proof is given in Appendix A.) Notice that, depending on a, jj assumes any value 
2 2 
between and fa. 

The consequence of Proposition 4 is that after an arbitrary rotation the appearance 
of an ellipsoid as it is approximated by the curvature method is in general neither better 
nor worse than its approximated appearance after a rotation about any of the main axes. 
As a result, if k models are required to cover all rotations around the X-axis, and k (or 
less) models cover all rotations around the X-axis, then at most k 2 models are required 
to cover all possible rotations in 3-D space. 


4 Summary 

The curvature method was designed to predict the appearance of objects with smooth 
bounding surfaces from different viewpoints. In this paper we have applied the method to 
ellipsoidal objects rotating arbitrarily in 3-D space and derived an expression describing 
the errors obtained, we have analyzed these errors and concluded that they depend on 
the shape of the ellipsoids, in other words, on the relative length of their axes, The 
error increases as the ellipsoid becomes “deep” (elongated in the Z-direction). We have 
shown that the errors are usually small, and that, in general, a small number of models 
is required to predict the appearance of an ellipsoid from all possible views. 

Finally, we would like to add that the curvature method discussed above is not re¬ 
stricted to contours originating from elliptic surface patches. It can as equally handle 
contours originating from hyperbolic patches - as long as the patches are visible. When, 
however, a patch is self occluded, a new aspect of the object is observed, and an addi¬ 
tional model should be utilized. The treatment of hyperbolic patches is demonstrated in 
Figure 6. Models of three tori with different radii were prepared analytically. The models 
were matched to an image that contained the tori in various positions and orientations 
using the curvature method. It can be seen that although the points of the inner circles 
of the tori come from hyperbolic patches, their prediction is still accurate. 
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Appendix A 


In this appendix we derive an expression of the error obtained when the curvature method 
is applied to a canonical ellipsoid rotating about the vertical axis. We then show that a 
similar error is obtained when the ellipsoid is rotating about any axis in space. Finally, 
we compute the error resulting from applying the curvature method to a non-canonical 
ellipsoid. 


A.l Rotation about the Vertical Axis 


Let 


x y z 

-h — H- 

a 2 b 2 c 2 


= 1 


be the surface of a canonical ellipsoid. Let pi = ( Xi,yi ) be a point on its silhouette. 
Assume the ellipsoid is rotating about the vertical (Y) axis by an angle 6. Let p 2 = 
(x 2 , 2 / 2 ) be the appeared position of pi following the rotation, and let p 2 = (ir 2 , y 2 ) be the 
approximated position of p 2 according to the curvature method. The relative error for 
the case of an ellipsoid that is rotating about the T-axis is given by 


E 


Proposition 3: The error is given by 


x 2 — x 2 
x\ 


E( — ,6) = cos 0 -\—-(1 — cos 0) — t/cos 2 9 H—- sin 2 9 


a 1 


a* 


a' 


Proof: The rim of a canonical ellipsoid contains the surface points for which z 

Therefore, the silhouette is defined by 

x 2 y 2 

+ ¥ = 1 

After the ellipsoid is rotated by an angle $ around the T-axis, it is described by 

(x cos 0 — z sin 0) 2 y 2 (x sin 9 + z cos 6) 2 


= 0. 


a‘ 


+ — + 
b 2 


= 1 


And its silhouette is given by 


x 


a 


2 cos 2 $ -f c 2 sin 2 9 b 2 


y 

+ — = 1 

•to 
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The position of p 2 = (:r 2 , 2 / 2 ) ls therefore 

x 2 = —\ja 2 cos 2 9 + c 2 sin 2 # 
a 

V 2 = Vi 


Next we calculate p 2 . Denote the surface of the canonical ellipsoid by the form 
F(x,y,z) = 1. According to Basri Sz Ullman [1988] the curvature radius with respect to 
the F-axis is given by 

F x c 2 x 

r* = -7 r = —r 

F zz a 2 

When the ellipsoid rotates around the F-axis by an angle 0, the position of p 2 is estimated 
by the curvature method to be 

x 2 — Xi cos 9 — —— 1 (1 — cos 0) 

a 2, 

h = y\ 

Consequently, the relative error is given by 

2 2 I 2 " * 

E( — ,6) = cos 6 H—-(1 — cos 6) — \ cos 2 6 H—- sin 2 9 

a 1 a 1 V a 2 


The error is therefore a function of 9 and 



A.2 Rotation in 3-D Space 

In this section we consider the case of a canonical ellipsoid rotating arbitrarily in 3-D 
space. A rotation in 3-D space can be decomposed into three successive rotations, about 
the Z-, Y- and Z-axes. The last rotation can be ignored since it does not deform the 
image and therefore does not change the errors. (The first rotation cannot be ignored 
since it determines the actual axis of the second rotation.) Let 

(—x sin a + y cos a 

b 2 

be the surface of a canonical ellipsoid rotated about the Z-axis by an angle a. We now 
examine this ellipsoid as it rotates about the F-axis by an angle 9. 

Proposition 4: The error is given by 

n2 

E <*•'> 



(a; cos a + y sin a)‘ 


a‘ 


+ 
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where 


C 2 
A 2 


c 2 c 2 

= — cos 2 a + — sin 2 a 
a 1 b 2 


Proof: In order to prove this proposition we have to show that every horizontal 

section of the ellipsoid defined above is an ellipse with an aspect ratio ^ as given in the 
proposition. 

Any nonempty intersection of an ellipsoid and a plane is either a point or an ellipse. 
The section is nonempty when y 2 < a 2 sin 2 a -f b 2 cos 2 a, and is a point when a strict 
equality holds. Given the canonical ellipsoid following its rotation about the Z -axis by 
an angle a, we show that the boundaries of its horizontal section can be represented as 

(g-Ip ) 2 z 2 

A 2 C 2 

which describes a canonical ellipse displaced along the X-axis. To establish the above 
relation, we show that for a constant value of y the surface equation of the rotated 
ellipsoid reduces to equation of the displaced ellipse. The two equations are identical if 
there exists a constant k ^ 0 such that the following equation system holds 


kC 2 
kC 2 x o 
kA 2 

kC 2 {A 2 - x 2 0 ) 


cos a 


+ 


sin 2 a 


a 4 


b 2 


,1 l \ 

y sin a cos £*( — - —) 


a 4 


1 


1-V( 


sin 2 a 


+ 


cos 2 a 


a 4 


b 2 


We obtain a system of four equations in four unknowns, A 2 , C 2 , x 0 and k. We now show 
that when y 2 < a 2 sin 2 a + 6 2 cos 2 a this system has a unique solution with positive values 
for A 2 and C 2 . 

Denote the right side of the four equations by 


P = 
<1 = 


r = 


s = 


cos 2 a 


sin 2 a 


a 2 + b 2 

.1 1 

y sin a cos ct{— --) 

b z a z 

1 




sin 2 a 


+ 


cos 2 a 


a A 


b 2 
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The solution to the system above is given by 


ar 0 = 


A 2 = 
C 2 = 
k = 


kC 2 x 0 q 

~kC 2 ~ = 7 
kC 2 (A 2 - x 2 ) + kC 2 x 2 

kC 2 

kC 2 (A 2 - x 2 ) 4- kC 2 x 2 
kA 2 

kA 2 p 2 r 
A 2 ps -f q 2 


ps + q 2 
P 2 

ps + q 2 
pr 


Notice that p,r > 0. This system therefore has a unique solution with positive values for 
A 2 and C 2 when ps + q 2 > 0. This inequality is satisfied when y 2 < a 2 sin 2 a A b 2 cos 2 a. 

Now, we can compute the value of the ratio ^ from this equation system by dividing 
the first equation by the third one 


C 2 c 2 2 c 2 2 
— = -cos a + -sin a 


Therefore, any horizontal section of this ellipsoid is an ellipse with an aspect ratio of 
§ 7 , and since translation does not affect the results of the curvature method, the error 
is given by 

fi2 

E (^*> 

Where ^4 2 and C 2 are the parameters of the ellipse, and 6 is the rotation angle around 
the Y- axis. 


A.3 Intermediate Models 


In this section we derive an expression of the error obtained when the curvature method is 
applied to an ellipsoid that is rotated about the T-axis (rather than a canonical ellipsoid). 
This computation is required for constructing Table 1 in Section 3.2. 

Let 

(a: sin a + z cos c*) 2 


(a: cos a — 0 sin a) 2 y 2 

r. I 1 r% I 


= 1 


a' 


b 2 c* 

be the surface of a canonical ellipsoid rotated about the T-axis by an angle a. Assume 
this ellipsoid is modeled by the curvature method. We consider now the error produced 
by using this model as the ellipsoid rotates about the T-axis by an angle 6. 
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Proposition 5: 


where 


The relative error is given by 
c 2 

E a (—,0) = cos 0 -f z' sin 6 + r'(l — cos 0) — x" 

(L 



Proof: Let pi = (x\,yi) be a point on the silhouette of the ellipsoid. Let Zi be its 

depth value, and let 7’j be its curvature value with respect to the K-axis. Then 

x 

X\ — 

a 

yi = y 

—x sin a cos a(a 2 — c 2 ) 

Zi = - 

a 

xac 2 

ri = - 3 

( a 2 cos 2 a + c 2 sin 2 0)2 




where p = ( x,y) is the corresponding point on the silhouette of the ellipsoid in its 
canonical position. 

Let p 2 = (a: 2 , 2 / 2 ) be the appeared position of pi after a rotation around the T-axis by 
an angle 0, p 2 is given by 

x 2 — —\Jcl 2 cos 2 (a -f 0) + c 2 sin 2 (a + 9) 

CL 

V2 = y 

Let p 2 = ( x 2 , y 2 ) be the position of p 2 approximated by the curvature method 

x 2 — X\ cos 0 Z\ sin 0 -)- (1 — cos 

fa = y 

Since yi = y 2 = y 2 , the error is defined by 

7-i x 2 - x 2 

H/qi — 

Xl 
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Let 



x 1 

£2 

Xl 


and we obtain the expressions given in the proposition. Notice that 

s 4' e) = E & e) 

in Eq. (2). 
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